On the falloff of radiated energy in black hole spacetimes 
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The goal of much research in relativity is to understand gravitational waves generated by a strong- 
field dynamical spacetime. Quantities of particular interest for many calculations are the Weyl scalar 
tpA, which is simply related to the flux of gravitational waves far from the source, and the flux of 
energy carried to distant observers, E. Conservation laws guarantee that, in asympotically flat 
spacetimes, 4'i oc 1/r and E ex 1/r^ as r — >■ oo. Most calculations extract these quantities at some 
finite extraction radius. An understanding of finite radius corrections to tp4 and E allows us to 
more accurately infer their asymptotic values from a computation. In this paper, we show that, 
if the final state of the system is a black hole, then the leading correction to 1JJ4 is 0{l/r^), and 



that to the energy fiux is 0{l/r 



not 0{l/r ) and 0{l/r ) as one might naively guess. Our 



argument only relies on the behavior of the curvature scalars for black hole spacetimes. Using black 
hole perturbation theory, we calculate the corrections to the leading falloff, showing that it is quite 
easy to correct for finite extraction radius effects. 

PACS numbers: 04.25. Nx, 04.30. Nk 
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I. INTRODUCTION 

Extracting radiation from the output of numerical cal- 
culations, as well as fluxes of quantities such as energy 
carried by radiation, is important for many problems in 
general relativity. Newman & Unti 1] provide an out- 
standing foundation for understanding analytically the 
asymptotic behavior of curvature tensors which deter- 
mine how gravitational radiation behaves as it propa- 
gates far from a radiating source. Perturbation theory 
also provides an excellent set of tools to help us under- 
stand the asymptotic behavior of radiation and fluxes. 

Many results on the distant behavior of radiation fields 
describe how quantities behave in the limit r — > 00. With 
the exception of characteristic methods (see, for exam- 
ple, [2[), most numerical calculations extract radiation at 
some large but finite radius r. Understanding the sub- 
leading corrections to the asymptotic behavior of radia- 
tive quantities could greatly improve our ability to ex- 
tract asymptotic fluxes and fields from numerical codes. 

Previous work ^^] found empirically that the form 



Eir) =E^{l + ^) 



does an outstanding job describing subleading corrections 
to the gravitational-wave energy flux. In this paper, we 
examine this behavior more carefully. In Sec. |lll follow- 
ing the formalism developed in Ref. [1|, we prove that 
this form is to be generically expected, and follows from 
the fact that at finite large radius r, the Weyl curva- 
ture scalar describing distant radiation takes the form 
i/,4(r) = -0°° (1 -h b2/r^). In Sec. Hill we use black hole 
perturbation theory to calculate the coefficients 62 and 
62. We conclude Sec. |IV]by discussing possible applica- 
tions of this result. 



II. TOOLS AND FORMALISM FOR 
UNDERSTANDING RADIATION FALLOFF 

A. Definitions 

We begin by defining the quantities which we will need 
for our analysis. Much of this discussion is adapted from 
Ref. [1|. We present these general definitions in some 
detail before specializing to the much simpler black hole 
case. 

Consider a vacuum, asymptotically flat spacetime. In- 
troduce a family of null hypersurfaces, each characterized 
by a constant parameter u. We take u = a;^ as one of the 
coordinates we will use to describe our geometry. Define 



L 



daU 



(2) 



Since these surfaces are null, the vector I" is tangent to 
null geodesies. This vector will be the first leg of a tetrad 
which we will use to characterize our geometry. Define 
r as the affine parameter along these geodesies; this will 
denote another of our coordinates. The remaining coordi- 
(1) nates x'^ {k € 3,4) then label the different null geodesies 
in each constant u hypersurface; they can be taken to be 
angles. 

We define a second null vector n" by requiring 



n'^ln 



1 



(3) 



To complete our tetrad, we next define a pair of unit 
spacelike vectors C" and p" that are orthogonal to ^", 
n", and each other. We then put 



(4) 
(5) 



m" = (C - zp")/V2 , 
m" = (C"+zp")/\/2 



We now use this tetrad to characterize the curvature 
of our spacetime. Let Ca^pu be the Weyl (vacuum) cur- 
vature tensor of the spacetime. Define the following 5 
complex Weyl projections: 



ip^ = —Cafi.pvn°''fh^n^'fff 



(6) 
(7) 
(8) 
(9) 
(10) 



Reference [1| shows that as we approach the asymptot- 
ically flat (r — )■ c») regime, these curvature components 
vary as follows: 



V-o = ^+0{l/r^), 



(11) 



^1 = ^+ ^'""^"^;:^"'''^°^ e^(iA^),(i2) 



V'3 



^3 f-a^As 



0{l/r^), 



(14) 



^. ^^-(^"^^^^;+^"^--^-)+e^(iA3).(i5) 

In Eqs. (IT2|) - ([T5|), the index fee [3, 4], the complex func- 
tion ^^ describes the angular components of the tetrad 
element to", and the functions ckRsc ^^nd 7r,sc are "Ricci 
spin coefficients," constructed by certain combinations 
and projections of the tetrad's covariant derivatives. For 
more details and discussion of these functions, see Refs. 
P,|j|. For our purposes, the most important fact to take 
from Eqs. (fTT|) - p5|) is that the leading falloff of ip^ is 
at 0{l/r). The subleading correction at ©(l/r^) is set 
by a coefficient that scales with A3, which controls the 
behavior of the curvature scalar -03. 



B. Perturbed black holes 

We now specialize to black holes. We use the Kinner- 
sley tetrad [^j which in Boyer-Lindquist coordinates is 
given by 



'" = ^(r2+a^l,0,a) , 



— (r2+a^~A,0,a 



\/2(r + iacos( 



(16) 

(17) 
(msin6',0, l,icsc6') . (18) 



For an unperturbed black hole spacetime, V'2 = — A//(r — 
iacos^)'^, and V'n = for n 7^ 2. Far from a perturbed 
black hole, '04 is also non-zero, describing the spacetime's 



outgoing gravitational waves: 

ipiir ^ go) = -[h+-ih^ 



1( 
2r V 



H. - iHy, 



(19) 



The Weyl scalar ipo is also generically non-zero for a per- 
turbed black hole, but we will not need its value in our 
analysis. Crucially for our argument, we can always put 
^3 = for our perturbed black hole Q. 

Comparing with Eqs. (fT2)) - (fTSi) , we read off 

(20) 
(21) 

(22) 



^3 = 0, 

A2 - -M 
1 



A4 



H, - iH, 



Combining these results with Eq. p^ . we see that cor- 
rections to ■04 come in at 0{l/r^), so that 



^.^^d^'^ 



(23) 



where 62 is a complex constant related to the (currently 
unknown) coefficient of this subleading falloff. 



C. Energy flux 

We now relate the curvature scalar 04 to the asymp- 
totic flux of radiation from the source. The energy flux 
in gravitational waves is given by 



E 



1 
16^ 



''dn 



{h+f + {h^f 



(24) 



Using Eq. ((T9)) . we can relate this to ^04 in the limit r 

cx): 
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i;°° = -^ lim 



'dfi 



dt ip4 



(25) 



Using Eq. ((23|) . let us now see what this implies about 
the behavior of E when radiation is extracted at some 
finite radius R. Let us first introduce a modal expansion, 
writing 



(26) 



For simplicity, we have taken the radiation to have a dis- 
crete frequency spectrum. The calculation can easily be 
extended to encompass a continuous spectrum. Combin- 
ing Eqs. dm and dH]), we find 



E{r) 



An ^-^ 



uj-'^ / r'^dn\iljX? 



'e:^ 



(27) 
(28) 



where 



S" = 



= {E-^r' X 



Attuj^ 



27rw2 



dn\A':, 



df7|A^|2(Re6^) 



(29) 



(30) 



In other words, an 0{l/r^) correction to V'4 produces an 
©(l/r"*) correction to E. We next must understand how 
to compute the coefhcient of this correction. We do so 
using black hole perturbation theory. 



III. SUBLEADING BEHAVIOR VIA 
PERTURBATION THEORY 

Perturbation theory is a powerful tool for calculating 
?/^4 and then determining fluxes such as E. In this sec- 
tion, we use black hole perturbation theory to confirm the 
general results of the preceding section, and to explicitly 
compute the magnitude of the subleading contributions 
to ip4 and E. 

Throughout this section we will assume a frequency- 
domain decomposition for ip^. This assumption means 
that solutions for ^^ separate Q: 



■4^4 



1 



imd) ^ — iujt 



(r — iacosi 



■^i?hn^(r)^,™(0;a^)e™^e 



The function Sim{9;auj) — S{9) is a spin- weighted 
spheroidal harmonic, and is discussed extensively in Ap- 
pendix A of Ref . [8] . It satisfies the eigenvalue relation 



10 de 



sm( 



^dS_ 
19 



(auj) cos" 9 + Aauj cos 
Am cos 9 + A 



sin^6l 



£ 



S^O. (32) 



In the a = limit, £—>/(? + 1), where / is the usual 
spherical harmonic index. 

The function Rimi^{r) = R{r) is governed by [3] 

A^|:(^^)-nOi?M--r(r), (33) 

often called the Teukolsky equation. Here and in what 
follows A = r^ — 2Mr + a} . Detailed discussion of the 
source Tir) is given in Refs. 0, Q- Fo'" our purpose, it 
suffices to note that an effective way to solve Eq. ([55)) is 
to first find a homogeneous solution, setting the source 
T(r) = 0. From these solutions, it is fairly simple to build 
a Green's function which we integrate over the source to 
find the particular solution for our problem. 

We show the potential V{r) in the Appendix. It 
depends on the eigenvalue £ via the parameter A — 
£ — 2amuj + o^lo^ — 2. An important property of V{r) 
is that it is long-ranged: as r — >■ cx), V{r) — )■ r^ . This 



makes computing R{r) for large r difficult. An excel- 
lent way to circumvent this difficulty is to first solve the 
Sasaki- Nakamura equation Q , 



^-F(O^-L.(r)X = 0, 

where r^{r) is the "tortoise coordinate," 
2Mr. , /r-r-\ 2i\/r_ 



(34) 



■In 



2M 



■In 



2M 



(35) 
The potentials F{r) and U{r) are also shown in the Ap- 
pendix. Their key property is that, unlike the Teukol- 
sky equation's V{r), they are short ranged: As r — ^ oo. 



F -^ (constant) /r and U 



(constant) /r"^. The 



solutions X{r) thus approach plane waves in the asym- 
potically fiat region, X{r -^ oo) — > gii^jr, _ Teukolsky 
equation solutions can be then be built from Sasaki- 
Nakamura equation solutions by the transformation 



R 



a - 



A 



X- 



Idx 
A dr 



(36) 



where x = XA/^/r"^ -\- a?. The functions a(r), /3(r), and 
7y(r) are listed in the appendix. 

A more accurate asymptotic form of X{r) is 



X(r) = v4°"*P°"'(r)e'"''- 
where 

[i/out 



A'"P'"(r)e- 



pin/out/ \ 



1 



Pi 



in/out 
P2 



in /out 

Pa_ 



(37) 



(38) 



LOT (ujr)^ (ujr)^ 

The coefficients appearing in this expansion are given by 
;,. i 



Pi 

pT 

and 



(A + 2 + 2amuj) , 



(39) 



- [(A + 2)2 - (A + 2)(2 - Aamcu) 
8 

4[arna; + SiMu — amLj{amLj + 2iMuj)\] , 



6 



[4amw + p'^i^ - 4 + 2amuj + 8iA/w) 



+12{Mujf - 2p'^\Muj - [aujf{\ - 3 
+2amu})] ; 



pT' = 


.,,.^ 


pT' - 


-^-^ 


pT ' 


= 4" + - 

Co 



rOJCi 



ijJ C2~ OJCi Pi + 



w^cs - ijj'^C2 (pi" -f- i) 
P2" + -^ - ^ + 2iMuj{aLom - 1) 



(40) 

(41) 

(42) 
(43) 



(44) 



The coefficients cq, Ci, C2, and C3 appear in tfie defini- 
tion of tlie function r](r), and are given in tfie Appendix; 
overbar denotes complex conjugate. 

The condition that radiation be purely outgoing far 
from the black hole picks out a solution of the form 



X = XooP°"*(r)e^' 



(45) 



as r — >■ 00. Performing the transformation (1361) . wc find 
that the Teukolsky solution R{r) can be written 

i?(r)==r3ZooQ™*(r)e'"^* (46) 

for r -^ cx), where Z^o — —'iXoo^'^ /cq, and where 



91 



Q°"*(r) = l + -i^ + 



92 



93 



(jjr (cjr)^ {ujrY 



(47) 



The coefficients 91,2,3 are given in the Appendix. 

Now use this solution to examine the fiux of energy a 
finite distance from the black hole. Using Eq. (l27l) . we 
find that 



Mr) = yj:r\Q°j'\' 



where 



\q: 



,out|2 



\z^ 



= 1 



47rLj2 ' 

6aaj(aa; — m) — X Af (A — 1) 



2u^r 



2^2 



(48) 

(49) 
(50) 



Notice that the leading correction to the energy fiux ap- 
pears at 0{l/r^), in agreement with Eq. (P5|. Compar- 
ing with Eq. ([T]), we find that the coefficient 62 which 
labels the 1/r^ falloff is 



62 



6auj{auj — m) — A 



(51) 



Recall that A is related to the spheroidal harmonic eigen- 
value £; ci. Eq. (15^ and following discussion. For 
Schwarzschild, this correction is particularly simple: 



62(0 = 0) = - 



l{l + l) 



2w2 



(52) 



where I is the spherical harmonic index associated with 
the mode under consideration. 



IV. DISCUSSION 

In this analysis, we have demonstrated that whenever 
one extracts radiation and radiative fluxes at a finite large 
radius, the subleading correction to these quantities is at 
an order ©(l/r^) beyond the leading asymptotic behav- 
ior. Hence, he correction to the curvature scalar tp^ is at 
©(l/r^), and to the energy fiux is at ©(l/r''). 

Using black hole perturbation theory, we have shown 
it is not difficult to calculate the coefficient of the sub- 
leading falloff, at least for a plane wave. The results we 



have found are consistent with the results shown in Ta- 
ble VI of Ref. Q. In that paper, a time-domain code 
was used to examine radiation from circular orbits. The 
time-domain code does not separate the angular behav- 
ior, and so many values of I are included in the analysis 
simultaneously. The radiation tends to be dominated 
by ^ = TO, with important but decreasing contributions 
from / = TO-|-1, I = m + 2, etc. Our expectation for the 
Schwarzschild radiation is thus likely to be close to the 
prediction from Eq. (j52p for I = m, skewed somewhat by 
contributions from I = m + 1. 

Let us test that prediction. Consider first the results 
for TO = 2. If we assume that the waves presented in Ref. 
3] for this case are dominated by radiation in the I = 2 
and / = 3 modes, then we expect 62 to be between 



62(0 = 0, / = TO, = 2) = 


-2..-^ , 


and (53) 


62(0 = 0,? = TO + 1 = 3) = 


-5uj-^ . 


(54) 



Table VI of Ref. [3] shows 

62(0 = 0, TO = 2) = -2.59uj~'^ 



(55) 



in reasonably good agreement with the intuition provided 
by our plane-wave expansion. Table VI also provides 
Schwarzschild data for m — 3; if those data are domi- 
nated hy I — 3 and I — 4, we expect 62 to be between 

62(0 = 0, / = r7i = 3) = —5lu~ , and (56) 

62(0 = 0,? = TO + 1 = 4) = -9cj"2 . (57) 

Table VI of Ref. [| shows 

e2(a = 0,TO = 2) = -6.20w~2 , (58) 

again agreeing reasonably well with the plane-wave ex- 
pansion. By computing the eigenvalues of the spheroidal 
harmonics for non-zero spin, one can likewise show that 
the Kerr values in Table VI agree reasonably well with 
the expectation of our plane-wave expansion. 

Bear in mind that the numerical magnitude of the cor- 
rection we derived strictly applies only for plane-wave 
expansions. As such, although we can provide good post 
facto justification of the coefficients of the subleading 
falloff, it would be difficult to predict those coefficients in 
advance. To do so, we would need to know the weighting 
of the different I modes which contribute to the radiation. 
Our only purpose in analyzing the coefficients shown in 
Ref. j^] is to show that the results presented there are 
consistent with our results here. 

For many calculations, it will not be worthwhile to de- 
compose the angular distribution of the waves, and thus 
to compute the subleading falloff in the manner shown 
here. It should be emphasized that the radial behavior 
of the falloff is independent of the I modes which con- 
tribute to the waves. As such, it would not be difficult to 
extract the radiation at several radii and simply fit the 
coefficient. That is what was done in Refs. Q and jlO| . 
Implementing such a multi-radius fit should make it pos- 
sible to more accurately extract the asymptotic radiation 



computed by numerical analysis, potentially reducing er- 
rors in such calculations by several percent. 

In general numerical spacetimes, it may be more com- 
plicated to take advantage of this result. The key ingre- 
dient to making the fallofF work as we have discussed is 
to choose a tetrad such that the Weyl scalar ^3 = 0. As 
long as one can perform a null rotation to put the space- 
time into such a "transverse" tetrad [y, \1^ , one should 
find find that subleading corrections to the flux of radia- 
tion fall off as 1/r^. It may be challenging to implement 
this rotation for the general case, but the improvement 
in accuracy could make it worthwhile. 
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and the function 

K{r) = (r^ + a'^)uj - ma . (A9) 

Recall that X — £ — 2amijj + a?uj^ — 2, where £ is the 
eigenvalue of the spheroidal harmonic. 

The potential V{r) appearing in the Teukolsky equa- 
tion (|33|) is given by 



Vir) 



K^ + 4i{r - M)K 



+ Siujr + A . (AlO) 



The potentials F(r) and U{r) appearing in the Sasaki- 
Nakamura equation (|M)) are 
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F{r) 



drj/dr A 
rj r'^ -\- a? 



(All) 



CM.#iM +M^-^„.)G„.) + c„f, 



(r^ + 0^)2 r^ + a^ 



where 



(A12) 



f/i(r) = V{r) + ^ 



dr \ A 



drj/dr ( dP / dr 
01 -\ 1 — 



(A13) 



Appendix A: Functions from black hole perturbation 
theory 

In this Appendix, we present various functions which 
arise in black hole perturbation theory that we need for 
our analysis. The functions ri{r), a(r), and /3(r) which 
appear in the transformation law (I36p are given by 



■n{r) 



Co + Ci/r + C2lr^ + c:i/r^ + Ci/r^ , (Al) 

4 + A . (A2) 
/3(r) = 2A [r - Af - 2A/r - iis:(r)] , (A3) 



tK{r)P{r) ^ ^.dK 
«M = A^ + 3*^ 



These functions in turn depend on the coefficients 

Co = -12?wM + A(A + 2) - l2aijj{aLo - m) , (A4) 

Ci = 8ia [iauj — X{auj — m)] , (A5) 

2 



C2 = ~2AiaM{auj - m) + 12a^ 1 - 2{auj - my 



C3 = 24ia^(aa; - m) - 2AM a^ , 
C4 = 120-^ , 



(A6) 
(A7) 
(A8) 



The coefficients p™2^ defined in Eqs. pop - (pi)) are 
found by requiring that the solution (1571) satisfy the 
Sasaki-Nakamura equation in each order in 1/r. After 
transforming to the Teukolsky equation solution R{r), 
the different orders in 1/r are labeled by the coefficients 
91,2,3 defined in Eq. (|T7| : 



?i 

92 



Pi 



93 



■™* i - ciuj/co , 



(A14) 



[-4c?a;2 -I- AcqU) [(p^* - i)ci + C20j] 



+cl [2 + 2ip7* - 4p^"' + A + Gamoj - UiMuj 
-6a V]] , (A15) 

-3 [4coCia;2 [(p™* - i)ci + 2c2UJ - 4c?a;=^] 
4cq 

+clu [^Alo [(p^' - i)c2 + Caw] -f ci [2 + 2ip°"* 

-4p^"* + A + Qamuj - 12iMuj - Ga^w^] ] 

+Co [4^3"' + 2t^ [M(5 + A) - bia^uj] 

-p^* [A - 2uj{3a^uj - Sam + AiM)] ] ] . (A16) 
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